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Abstract 

Wc investigate a lattice model of polymers where the nearest-neighbour monomer-monomer 
interaction strengths differ according to whether the local configurations have so-called 
"hydrogen-like" formations or not. If the interaction strengths are all the same then the 
classical 0-point collapse transition occurs on lowering the temperature, and the polymer 
enters the isotropic liquid-drop phase known as the collapsed globule. On the other hand, 
strongly favouring the hydrogen-like interactions give rise to an anisotropic folded (solid- 
like) phase on lowering the temperature. We use Monte Carlo simulations up to a length of 
256 to map out the phase diagram in the plane of parameters and determine the order of 
the associated phase transitions. We discuss the connections to semi- flexible polymers and 
other polymer models. Importantly, we demonstrate that for a range of energy parameters 
two phase transitions occur on lowering the temperature, the second being a transition from 
the globule state to the crystal state. We argue from our data that this globule-to-crystal 
transition is continuous in two dimensions in accord with field-theory arguments concerning 
Hamiltonian walks, but is first order in three dimensions. 
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1 Introduction 



The self- avoiding walk (SAW) on a lattice [U [2] is a key model in statistical mechanics for the 
study of the static properties of polymers. Incorporating interactions in this model make it 
possible to represent many features of real polymers. Regardless of the constraints of the lattice 
relative to the real world, the model mimics very well many properties of physical systems [3]. 
The self-avoiding walk on a lattice is a random walk which is not allowed to visit a lattice site 
more then once. Each visited lattice site is considered as a monomer of the polymer chain. A 
common way [2] to model intra-polymer interactions in such a walk is to assign an energy to each 
non-consecutive pair of monomers lying on the neighbouring lattice sites. This is the canonical 
ISAW model which is the standard model of polymer collapse using self-avoiding walks. With 
this modification one studies a polymer in a solvent, where the energy between monomers can 
be attractive or repulsive and depends on temperature. If the energy is repulsive the polymer 
behaves as a swollen chain (the so-called excluded-volume state) regardless of temperature and 
one says that it is in a good solvent. When the energy is attractive, and the temperature is 
low enough, the chain becomes a rather more compact globule 0], reminiscent of a liquid 
droplet: this is also known as the poor solvent situation. The transition point between those 
two phases is called the #-point; it is a well studied continuous phase transition (see [5] and 
references therein). 

The modelling changes as soon as we want to describe any biological system (e.g. proteins), in 
which the hydrogen bonding plays an important role [B|. One of the main features of the bonding 
is that the interacting residua lie on a partially straight segments of the chain. Hydrogen-like 
bonding was first modelled on the cubic and square lattices using Hamiltonian paths by Bascle et 
al. [7J. A monomer acquires a hydrogen- like bond with its (non-consecutive) nearest neighbour 
if both of them lie on straight sections of the chain (see Figure [2]) . The interacting self-avoiding 
walk modified to have only such interactions will be referred to as the hydrogen-like bonding 
model, or rather hb-model. The hb-model was studied in mean-field approximation [7J and a 
first-order transition from a high-temperature excluded-volume (swollen) phase to a quasi-frozen 
solid-like phase was found in both two and three dimension. Hence this would indicate that it 
is a different type of transition from the #-point. Note also that the low temperature hb-phase 
is anisotropic whereas the collapsed globule of the standard #-point model is isotropic. The 
hb-model on the square lattice was studied directly by Foster and Seno by means of the transfer 
matrix method [8] and by Krawczyk et al. [9] on both the square and cubic lattice using a Monte 
Carlo method. In both studies a first-order transition was found between an excluded-volume 
(swollen-coil) state and an anisotropic ordered compact phase in two and in three dimensions, 
again in opposition to the #-point [2]. 

It is appropriate to compare this difference between the behaviour of the hb-model and 
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#-point models with the difference between the behaviour of interacting semi-stiff polymers 
and the fully-flexible 0-point polymers. This is because hydrogen bonding induces an effective 
stiffness in the polymer between those monomers that are taking part in the interactions. As the 
temperature is lowered the proportion of the monomers experiencing this stiffness increases so 
while not all the segments of the polymer feel this stiffness at high temperatures, the proportion 
of monomers involved with nearest-neighbour /16-interactions increases towards unity as the 
temperature is lowered. In three dimensions, Bastolla and Grassberger [10] discussed so-called 
semi-stiff self-avoiding walks, which interact via all nearest-neighbours, as in the (9-point model, 
and include a bending energy. They showed that when there is a strong energetic preference for 
straight segments, this model undergoes a single first-order transition from the excluded-volume 
high-temperature state to a state similar to the low-temperature solid-like state of the hb-model. 
Intriguingly, if there is only a weak preference for straight segments, the polymer undergoes two 
phase transitions: on lowering the temperature the polymers undergoes the #-point transition to 
the liquid globule followed at a lower temperature by a first-order transition to the frozen phase. 
We should point out though that in two dimensions the transition between the globule and the 
frozen state has only been studied in Hamiltonian walks, and there it seems to be continuous 
one [TT] . 

To complicate matters further, there is at least one other model using a different definition 
of interactions which could be regarded as hydrogen-like bonding. This model [12] defines 
interactions between parallel segments, that is, bonds of the lattice occupied by the walk and 
so connecting monomers, see Figure [2j We will call this model the interacting bond model. 
Studying this model by means of Bethe approximation, Buzano and Pretti [12] found, in both 
two and three dimensions, two phase transitions: while decreasing the temperature the ^-collapse 
to isotropic globule phase is followed by a first order transition to a solid-like phase. Hence, 
this is similar to the semi-stiff model for weak stiffness. The interacting-bond model in two 
dimensions has recently been studied by Foster [H] and also displays two transitions. In [13] 
Buzano and Pretti added isotropic nearest-neighbour monomer-monomer interactions to the 
interacting-bond model, and investigated the phase diagram in three dimensions, again in the 
Bethe approximation. They showed that the phase diagram is similar to the interacting-bond 
model. However, if the interaction between monomers are repulsive, there is only one first-order 
phase transition from the swollen coil to the solid-like phase. This is again reminiscent of the 
semi-stiff model for strong stiffness. 

It is therefore of some interest to study an enhanced hb-model where non-hydrogen bond near- 
est neighbours are also considered. Hence, in this paper we investigate a model of self-interacting 
self-avoiding walks with two types of nearest-neighbour interaction: the /^-interactions and 
nearest-neighbour interactions that are not hydrogen bonds, which we denote as n/i-interactions. 
The competition between these two types of interaction (hb versus nh) leads to a three-phase 
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phase diagram, with excluded volume, globule and frozen phases. Of special interest is the com- 
parison to the semi-stiff model. One key question is whether there can exist two phase transitions 
on lowering the temperature. The order of the transitions in two and three dimensions are also 
of interest. We use a Monte Carlo technique, known as FlatPERM [15], to study self-avoiding 
walks on the simple cubic and square lattices with interactions as described. 

The paper is organised as follows. In Section 2 we explain more carefully details of the 
model. In Section 3 the phase diagram in both dimensions is discussed. A discussion of the 
anisotropy of the model is also given. We conclude with a summary and discuss the similarity 
of this model to the semi-stiff model. 

2 Model and simulations 

The polymer is modelled on a square and simple cubic lattice as a self-avoiding walk with 
interactions between different types of nearest-neighbours monomers: that is monomers that 
are not consecutive in the walk but nearest neighbours on the lattice. The strength of the 
interaction depends on the relative position of the monomers involved in the interaction to 
those next to them on the walk. A segment is defined as a site along with the two adjoining 
bonds visited by the walk, and we say that a segment is straight if these two bonds are aligned. 
The 'hydrogen bonds' are nearest neighbour interactions that are between monomers where both 
are part of straight segments of the polymer, see Figure HJ 




Figure 1: The types of nearest-neighbour interactions between two straight segments of the 
polymer involved in the /^-interactions: parallel segments (left) and orthogonal segments (right). 
In the model studied in this paper these two types are weighted equally. In two dimensions, 
only parallel interactions are possible. 

Our model weights both the parallel and orthogonal hb interactions equally. Our model also 
includes all other possible nearest neighbour interactions and assigns them a different Boltzmann 
weight. The two types of distinguished interactions are shown in Figure [2j As just described 
interactions between monomers sitting on straight lines, as in monomer #5 of the triple of 
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monomers #2, #5 and #14, form the hydrogen-bond interactions with like monomers, as in 
the interacting pairs of monomers 2-5 and 5-14. The energy of those interactions is denoted as 
—£hb- The other kind of interaction consist of any nearest-neighbour interaction between non- 
consecutive monomers that are not hydrogen bonds, as in the pairs 1-6 and 10-13 for example. 
The energy of those interactions is denoted as — £nh- 
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Figure 2: Interactions in 3d interacting self-avoiding walk. In our model we distinguish two differ- 
ent kind of interactions, which are denoted by two different colours. The red colour (dark shad- 
ing) denotes hydrogen-like interactions (hb) (monomers 2-5, 5-14) whereas the green colour (light 
shading) denotes all other interactions between two neighbouring monomers (nh) (monomers 
1-6 and 10-13 for example). 

Note that the number of all nearest-neighbours interactions, m, is equal to the sum of the 
number of the two types of interaction considered in our model, that is m = m n i + m^. The 
energy of configuration (p n of an n-step walk is calculated as 

E((p n ) = -m h b((p n ) ■ e hb - m nh (ip n ) ■ e nh , (2.1) 

where rrihb and m n h are the number of hydrogen-like bond interactions and non hydrogen-like 
nearest-neighbour interactions, respectively. The inverse temperature is denoted as (3 = 1/ksT, 
where ks is the Boltzmann constant and T the absolute temperature. We define for convenience 
Phb = @£hb and j3 n h — 0£nh- The partition function is then given by 

Z n (Phb,Pnh)= E C n,m hb ,m nh e^ m "+^~ (2.2) 
rn hb ,m nh 

with C njmhb)rnnh the density of states. Canonical averages are calculated with respect to this 
density of states. 

Since we will consider simulational results along lines in parameter space at a constant ratio 
of e nh /£hb = Pnh/Phb we define e nh = 7, and e hb = 1 - 7, so that the energy is then given by 

E = -m hh • (1 - 7) - m nh • 7. (2.3) 
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In our study we will analyse the (reduced) specific heat to investigate the phase diagram: that 
is, 

_ 1 (BVW ( , 4) 
1 n 

We will also consider fixing one of the parameters either flhb or (5 n h and varying the other. 
To analyse the possible phase transitions we then use the fluctuations in number of monomers 
of appropriate type. In the case of /3hb being constant we consider 

o 2 {m n h) = (ml h ) - (m nh ) 2 . (2.5) 

When f3 n h is constant we consider 

v 2 {m h b) = (m 2 hb ) - (m hb ) 2 . (2.6) 

Simulations have been performed with FlatPERM algorithm [15j . We have simulated the 
models using a two parameter implementation (utilising mhb and m n h) for length n = 128 where 
the simulation directly estimates this density of states C ntTnhbtmnh . We have also performed one 
parameter {mhb or wink) simulations for systems of size 256 where the simulation estimates 
partial summations of this density of states over one of the variables. 



3 Results and discussion 

3.1 Pure hydrogen bonding and the canonical IS AW models 

When 7=1/2 the model becomes the canonical interacting self avoiding walk (IS AW), which 
displays the ^-transition from coil to globule state. The ^-transition is a second-order phase 
transition in both two and three dimensions. In two dimensions the established crossover expo- 
nent is <p = 3/7 [16] . which implies a negative specific heat exponent a = 2 — l/(p = —1/3 (the 
specific heat does not diverge on approaching the transition) [T7]. In three dimensions, which is 
the upper critical dimension for the ^-transition, the specific heat is expected to diverge loga- 
rithmically [18] . The collapsed state is an isotropic dense liquid- like droplet with a well-defined 
surface tension [191 [20] . 

On the other hand, for 7 = the model becomes the hb-model studied by Foster and Seno 
[S] on the square lattice and Krawczyk et al. on the square and simple cubic lattices [U] • In both 
two and three dimensions, there is a single first-order transition to a folded crystalline state, 
which is anisotropic. 

Using these results as a starting point, one therefore expects there to be at least these three 
phases (swollen, globule, crystal) in the full two-dimensional parameter space explored here. 
For 7 > 1/2, the /ifr-interactions are suppressed relative to the non-hb interactions, and the 
simplest hypothesis would be that the ^-transition is not affected. To test this hypothesis, we 
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have considered the line 7 = 1 below. For 7 close to zero, the /^-interactions dominate, and one 
may expect that there exists some range of values for 7 for which the first-order transition of 
the pure hb-model persists. To test this, we need to consider a small value of 7. We of course 
then need to consider other values of 7 to see if the two transitions can occur for fixed 7 and 
whether there exist any other phases. 

3.2 Features of the phase diagram 

To gain an understanding of which values of 7 we may need to consider more closely, we first 
examine the fluctuations in the numbers of interaction across a wide range of (PhbiPnh)- As 
in previous work [21|, [2~2t [9] , we found the use of the largest eigenvalue of the matrix of second 
derivatives of the free energy with respect to the parameters f3 n h and f3hb most advantageous 
to show the fluctuations in a unified manner. Figure [3] displays density plots of the size of 
fluctuations for —0.1 < (3 n h,Phb < 2.0 in two and three dimensions, respectively. The lighter the 
shade the larger the fluctuations. 

It suggests the presence of three thermodynamic phases separated by three phase transition 
lines meeting at a single point. From our discussion above, we can therefore identify these three 
phases as swollen, globule, and crystal. Therefore, for small values of fi n h and Phb, we expect the 
model to be in the excluded volume universality class of swollen polymers, since at /3 n h = Phb = 
the model reduces to simple self-avoiding walks. For fixed (3hb an d large (5 n h, we deduce that the 
polymer is in a globular state, while for fixed (3 n h and large (3^, we deduce that the polymer is 
in the anisotropic crystalline state. 

A fixed value of 7 corresponds to a straight line of slope (1 — 7V7 on the plots in Figure El 
The line with 7 = 0.5 is shown on both plots and it is clear that for any value of 7 > 0.5 (slope 
less than 1) the system will undergo only one phase transition on lowering the temperature. 
This transition should be in the universality class of the 0-transition. Below we consider the 
line 7 = 1 to verify this. The line with 7 = 0.1 is shown on both plots and it is clear that 
there is indeed a range of values of 7 around (slope sufficiently large) for which the system will 
undergo a single first-order hb-model like transition on lowering the temperature. The figures 
also suggest that there exists a critical value of 7, say, j c , where this scenario ends. 

Therefore one may deduce that there exists some range of 7 between 7 C and 0.5, not neces- 
sarily the whole range, for which the system undergoes two transitions on lowering the temper- 
ature. In this region, the polymer starts in the swollen state at high temperatures, undergoes 
a ^-transition to a globular state on lowering the temperature, and on lowering the tempera- 
ture further, undergoes a further (novel) transition to the crystalline state. We verify this by 
considering the line with 7 = 0.4. 

Figure shows the specific heat in two and three dimension for n = 128 at 7 = 0.1,0.4, and 
1.0 as a function of temperature. For 7 = 0.1 and 7 = 1.0, there is one peak in the specific heat 
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Figure 3: Density plots of the logarithm of the largest eigenvalue of the matrix of second deriva- 
tives of the free energy with respect to (5 n h and 0hb (the lighter the shade the larger the value). 
Both plots are for n = 128, for two- and three-dimensional systems, top and bottom, respectively. 
The lines shown indicate cross-sections for which we have performed additional simulations or 
analysis. In both pictures we show lines with slope (1 — 7V7 for 7 = 0.0,0.1,0.4,0.5 and 1.0. 
Displayed are also vertical lines at (5 n h = 1.0 and 0.7 in two and three dimensions, respectively. 

which is sharp for 7 = 0.1 and relatively broad for 7 = 1.0. This is consistent with the scenario 
described above, where at 7 = 0.1 there should be a first-order transition in the thermodynamic 
limit, while at 7 = 1.0 we expect a transition in the #-point universality class. Also as predicted 
above, for 7 = 0.4 there are two well-formed peaks in the specific heat. In three dimensions the 
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Figure 4: The specific heat for different values of 7 in two and three dimensions for n = 128. 
For 7 = 0.1 and 1.0 we see only one maximum in the specific heat, so we expect only one phase 
transition. For 7 = 0.4 two maxima occur: upon changing the temperature there are two phase 
transitions. (The specific heat for 7 = 0.1 is divided by a factor 15 in three dimensions and 10 
in two dimensions to allow to depict all curves on one diagram.) 

peak at lower temperature is sharp, while the one at higher temperatures is relatively broad, 
consistent with a #-like transition from the swollen coil to a collapsed globule at a moderate 
temperature, followed by a stronger globule-crystal transition at a lower temperature. In two 
dimensions, there are two peaks of roughly equal height. However, they are not well separated, 
which indicates that we need to go to longer lengths to study these transitions. 
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3.3 Low temperature phases 

Before considering the order of the phase transitions, especially the globule-to-crystal transition, 
we verify that the low-temperature phases have the properties assumed above. In particular, we 
demonstrate that while the globular phase displays no orientational order, the phase for large 
/3hb at fixed (5 n h is a crystal phase which displays strong orientational order by showing that in 
this phase the bonds between monomers prefer to align with one axis of the lattice. 
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Figure 5: The plots display the anisotropy parameter p of the system in two and three dimen- 
sions. 



To detect orientational order, we utilize an anisotropy parameter defined in Bastolla and 
Grassberger [10]. If we denote the number of bonds parallel to the x-, y-, and z-axes by n x , n y , 
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and n z , respectively, we define 



p = 1.0- 



mm(n x ,n y ,n z ) 
max(n x ,n y ,n z ) 



(3.7) 



In a system without orientational order, this quantity tends to zero as the system size increases. 
A non-zero limiting value less than one of this quantity indicates weak orientational order with 
n min oc n max , while a limiting value of one indicates strong orientational order, where n max ^> 



We consider a fixed value of (3 n h such that the system is collapsed for any value of (3hb- For 
small values of fihb the polymer is in the globular phase, while for large values it is expected 
to be in the crystal phase, see Figure El In two dimensions, we use i3 n h = 1.0, while in three 
dimensions we use (3 n h = 0.7. 

Figure [5] shows p as a function of Phb for different lengths ranging from 32 to 128 in two 
and three dimensions. For small 0^, we find that p converges to zero as n -1 / 2 as expected 
if only statistical fluctuations are present. Similarly, for large (3^, we find that p converges 
to one in a corresponding fashion. This indicates the presence of strong orientational order 
in the large /3/jfe-phase, which we then deduce to be the ordered crystal. Intriguingly, for two 
dimensions only, there exists a value of Phb, 1-54, at which p seems to be independent of system 
size. Moreover, by choosing an appropriate exponent {(f) ~ 0.75), one can show indications of a 
scaling collapse of the data near this point. This is our first indication that the two-dimensional 
and three-dimensional globule-to-crystal transitions are different. 

3.4 The phase transitions 

3.4.1 Swollen coil to folded crystal and swollen coil to globule 

As discussed elsewhere [9], when 7 = 0, that is (5 n h = 0, there is a single phase transition which 
is first-order in both two and three dimensions. We have verified that for small 7, for which we 
choose 7 = 0.1, this scenario remains intact. A bimodal distribution can be clearly seen forming 
as the system size becomes larger very strongly in three dimensions, and more weakly in two 
dimensions. 

To show that the ^-transition extends from 7 = 1/2 to larger values of 7, we can focus on 
the case of 7 = 1, which means that P^b = 0.0 and /^-interactions are irrelevant. 

In two dimensions the maximum of the fluctuations per monomer for length n for the #-point 
behaves as 



where (j> = 3/7. We have estimated from our data collected at short length that <j) = 0.49. This 
is consistent with the observation that for finite system size the effective crossover exponent 
decreases from a value well above 0.5 [5] to the theoretical value predicted for the #-point 
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2d 




cj) = 3/7. Having extrapolated the limiting value ^(m,^), we show in Figure [6] a scaling plot 
of the dependence of the singular part of the fluctuations as a function of f3 n h- 

In three dimensions [2j [18] theory predicts 4> = 1/2 and a logarithmic divergence of the 
maximum of fluctuations. The scaling of the fluctuations around the transition for P^b = are 
shown in Figure El At short lengths, we find strong corrections to scaling, in accord with the 
observations in [23\ 124] . The effective exponent of the logarithm is equal to 2.8, which is about 
an order of magnitude more than the value 3/11 predicted [18j . However we have checked that 
the exponent decreases with the system size. 
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3.4.2 Collapsed globule to folded crystal 

We begin with the three-dimensional case and return to our simulations at fixed j3 n h = 0.7. 
On varying (3hb we find a first-order transition from the globule to the crystal. The maximum 
of fluctuations per monomer cr 2 (m/ l {,)/n increases linearly in n, and the shift of the inverse 
temperature scales as 1/n, in accord with finite-size scaling of a first-order transition. Figure [7] 
shows the corresponding scaling collapse, with an extrapolated value of /3£ b = 1.34. 

We turn to our simulations of the two-dimensional case at fixed (3 n h = 1.0. Now, on varying 
/3hb we find a transition which is much stronger than the #-point, but shows no indication of being 
first-order: the maximum of fluctuations per monomer <T 2 (m/ l fe)/n diverges with an exponent 
less than one. A scaling plot using a consistent power law is not convincing. Our best estimate 
for the crossover exponent comes in fact from the scaling of the anisotropy parameter p discussed 
above, which gives (3^ = 1.54 and a crossover exponent in the vicinity of <p = 0.75. The scaling 
of the fluctuations is not inconsistent with these values. 
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Figure 7: A scaling plot for the fluctuations in rrihb for the globule-to-crystal transition occuring 
along the line (3 n h = 0.7 in three dimensions. 



3.5 Summary 

We summarize our findings by presenting conjectured phase diagrams in two and three dimen- 
sions in Figure [8l For large values of the ratio of the interaction strength of hydrogen-bonds to 
non-hydrogen bonds, a polymer will undergo a single first-order phase transition from a swollen 
coil at high temperatures to a folded crystalline state at low temperatures. On the other hand, 
for any ratio of these interaction energies less than or equal to one there is a single #-like tran- 



13 



2d 



crystal 




swollen 
coil 




compact 
globule 






Pnh 






3d 


crystal 




i 
i 

swollen I 
coil i 




compact 
globule 


Pnh 



Figure 8: The schematic conjectured phase diagrams in two and three dimensions. A solid 
line denotes a first-order transition, and a dashed line a #-like transition. The dot-dashed line 
represents a putative second-order phase transition that is not yet completely characterised. 

sition from a swollen coil to a liquid droplet-like globular phase. For intermediate ratios two 
transitions can occur, so that the polymer first undergoes a 0-like transition on lowering the 
temperature, followed by a second transition to the crystalline state. In three dimensions we 
find that this transition is first order, while in two dimensions we find that this transition is 
probably second order with a divergent specific heat. It can be argued using a zero-temperature 
argument that this scenario exists as soon as the ratio of the interaction energies is greater than 
one. In this way the phase diagram described is qualitatively similar to the one of the semi-stiff 
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interacting polymer model described in three dimensions by Bastolla and Grassberger |10j . 

The interesting questions that remain for future work include further characterizing the 
globule-to-crystal transition in two dimensions and also clearly delineating the range of interac- 
tion ratios for which two transitions occur. 
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